IMPROVED MAPPING OF cc ONTO CHARMONIUM STATES * 
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We discuss the relative rates of different charmonium states and introduce an 
improved model for mapping the continuous cc mass spectrum on the physical 
charmonium resonances. 



The theoretical description of charmonium production separates the hard 
and soft parts of the process based on the factorisation theorem in QCD. 
Perturbative QCD (pQCD) should be applicable for cc production, since the 
charm quark mass m c is large enough to make a s (m^.) a small expansion pa- 
rameter. The formation of bound hadron states occurs through processes with 
small momentum transfers making a s large and preventing the use of pertur- 
bation theory. This forces us to use phenomenological models to describe the 
formation of charmonium states from perturbatively produced pep. pairs. The 
Colour Evaporation (CEM) 111 and Soft Colour Interaction (SCI) &B models arc 
based on a similar phenomenological approach, where soft colour interactions 
can change the colour state of a cc pair from octet to singlet. They employ the 
same hard pQCD processes to produce a cc pair regardless of its spin state. 
A colour singlet cc pair with an invariant mass below the threshold for open 
charm (m cE < 2m£>) will then form a charmonium state. Mapping cc pairs 
onto charmonium states is based on phenomenological parameters taken from 
comparison to data (CEM), or from spin statistics (SCI), the latter resulting 
in a fraction of a specific quarkonium state i with total angular momentum Jj 
and main quantum number rij given by ft = , where T = (2 J, + l)/rij. 

The CEM and SCI models give in general a good description of charmonium 
data both from fixed target and Tevatron energies BaB. However, the over- 
all energy-independent factors which describe the mapping of cc pairs onto 
different charmonium states are not satisfactorily understood. Jn particular, 
spin statistics gives an acceptable ratio ip'/tp at the Tevatron 0, but fails at 
fixed target energies where an extra suppression by a factor four is requireda. 
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In other words, this ratio has an energy dependence which has not been ac- 
counted for. In this work we try to understand this problem by developing 
a model which introduces a correlation between the invariant mass of the cc 
pair and the mass of the final charmonium state. 

The cc pair is produced in a pQCD process with a continuous distribution 
of its invariant mass m c g and must be mapped onto the discrete spectrum of 
charmonium states. The soft interactions that turn the pair into a colour sin- 
glet may very well change its mass by a few hundred MeV, which is the typical 
scale of the soft interactions, but larger mass shifts should be suppressed. The 
different charmonium states (?y c , J/tp, Xcj ^' ) are i however, separated in mass 
over a region of about 1 GeV. It therefore seems likely that the probability to 
form a particular charmonium state will depend on the original value of m cS 
and not just on an overall factor, such as spin statistics. Thus, there should 
be a relatively larger probability to form states that are nearby in mass, than 
those further away on the mass scale. Based on these considerations we have 
constructed the following model. 

The smearing of the cc mass due to soft interactions is considered as 

n i \ ( (mcz-m) 2 \ 

Gs m e(m c -,m) = cxp — — I, (1) 

V ^ G sme J 

where the gaussian width a sme should be a few hundred MeV. The different 
charmonium resonances i are represented by functions Fi(rrii, m) around their 
physical masses m*. For a cc pair with a given mass m cS , the probability to 
form a certain charmonium state is then 

p ^ m s_ J G sm e(rn C c,m)Fi(mi,m)dm ^ s,G sm eKc, ^) 

m)Fj(mj,m)dm £\ SjG sme (m c -, rrij) ' 

where the sum is over all charmonium resonances. The approximate result, 
which is used in our practical calculations, is obtained by letting Fi(m,i, m) = 
Si5(m — nii), i.e. neglecting the width of the very narrow charmonium states 
and including the previously used spin statistical factor Sj = (2Jj + l)/rij. 

The smearing of m cE across the threshold 2mo for open charm, implies 
non-zero contributions for charmonium also above the DD threshold as well 
as some open charm production for m cE originally below this threshold. For 
a given original m cg value, the probability to end up in the region above DD 
threshold is given by the area of the smearing gaussian in that region, i. e. 

A{m C c) = — == / exp — dm . (3) 



2na sm e J2m D V 

The probability to end up in a specific resonance is then given by (1 — A)Vi 
and is shown in Fig. [l] for the different charmonium states. For a given original 
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Figure 1. Probability distributions for the different charmonium states as obtained in the 
model with gaussian smearing (<r sme = 400 MeV). The resulting total probability for 
charmonium production and the remainder as open charm production are also shown. 

cc mass m cE one can here see the fractional production of different states, as 
well as the sum of all charmonium states together with the remainder giving 
open charm. 

By folding these charmonium probability functions with the distribution 
in m cg obtained from pQCD, i.e. da C c/dm C c (LO processes and m c = 1.35 
GeV) , the cross section for a given charmonium state is obtained as 



Applying this mapping procedure to the CEM model we obtain the results 
in Fig. H As opposed to the simple spin statistics factor, this model gives a 
reasonable description of the observed ratio of ip' to J / ip production and frac- 
tions of J/ip produced directly, coming from decays of \c states and from ip'. 
In particular, the model gives a characteristic energy dependence of the kind 
indicated by the data. Considering all observables from Fig. ^, the preferred 
value for the a sme is 400 MeV or lower. This fits well with the expectation for 
the non-perturbative dynamics that this model should describe. The essen- 
tial new ingredient of this model, namely introducing a correlation between 
the invariant mass of the cc pair and the masses of the different charmonium 
states, is found to give an improved description of the ratios of different char- 
monium states as compared to using constant factors based on spin statistics. 
In particular, these ratios acquire an energy dependence which is particularly 
strong at low energies where threshold effects are more pronounced. 
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Figure 2. The ratio ip' / J /if) (times their branching ratios for decay into fi + fi~) (left) and 
fractions of J/i/> produced directly, and coming from the decay of Xc and ip' states (right) in 
hadron-hadron interactions of cms energy i/s. Data 6 ' 7 - 8 compared to simple spin statistics 
and to our model with different gaussian smearing widths applied to CEM. 
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